function BootStatsIS = BootIS(XBoot,YBoot)

global BootISstruct;

[T, NBoot] = size(YBoot); % XBoot is assumed to be only a single variable
BootStatsIS = repmat(BootISstruct,NBoot,1);
for b = 1:NBoot
    y = YBoot(2:end,b);
    Xb = XBoot(:,b); Xb = (Xb - mean(Xb))/std(Xb);
    x = [ones(T-1,1) Xb(1:end-1)];
    [beta,tOls,~,tNwest,~,adjr2] = ols(x,y,1);
    BootStatsIS(b).beta = beta(2);
    BootStatsIS(b).tOls = tOls(2);
    BootStatsIS(b).tNwest = tNwest(2);
    BootStatsIS(b).adjr2 = adjr2;
end
return